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Abstract 

The stability of (/-Gaussian distributions as particular solutions of the linear diffusion equation 
and its generalized nonlinear form, ^^g^'*^ = ^^^dJ^^ — porous-medium equation, is inves- 
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in such a way that the relaxation rule for the kurtosis evolves in time according to a g-exponential, 
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an infinite- variance distribution (q^ > 5/3) into a finite- variance one {q < 5/3), the relaxation 
towards the asymptotic solution may occur very slowly in time. This fact might shed some light 
on the slow relaxation, for some long-range-interacting many-body Hamiltonian systems, from 
long-standing quasi-stationary states to the ultimate thermal equilibrium state. 
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I. INTRODUCTION 



The linear diffusion equation, which is one of the most important differential equations 
of classical physics, rules the time evolution of the probability distribution associated with a 
particle diffusing in a homogeneous medium. This probability distribution spreads in time, 
with its second moment increasing linearly with time, being appropriate for a description of 
many physical phenomena, usually classified as normal diffusion. Such an equation presents 
a structure very common in physics, e.g., it is analogous to the heat-conduction equation, 
and by introducing an additional term, characterizing an external harmonic force field, it 



becomes the linear Fokker-Planck equation (FPE) [l], O]. The linear FPE is essentially 
associated with the Boltzmann-Gibbs (BG) formalism, in the sense that the Boltzmann 
distribution, which is usually obtained through the maximization of the BG entropy under 
certain constraints, also appears as the stationary solution of the linear FPE . However, 
such simple linear equations are not appropriate for dealing with many physical situations, 
characterized by anomalous diffusion, like particle transport in disordered media 4] and 
motion in optical lattices j^. 

Due to the recent advance in computer technology, some problems in physics that re- 
mained unexplored for a long time are now under investigation, at least from the computa- 
tional point of view. In particular, one should single out those problems described in terms of 
nonlinear differential equations, which have led to new features and interesting puzzles that 
keep challenging many physicists nowadays. The nonlinear FPEs {g], which are formulated 
by introducing modifications in the standard FPE, appear naturally as good candidates for 
describing anomalous-transport processes. In most cases, the nonlinear FPEs are proposed 
as simple phenomenological generalizations of the linear FPE 
sible to obtain them through approximations in the master equation 



|, alt houg h it is pos- 
n, 3, 13|. Recently, 



a relation involving quantities of the FPE and entropic forms was proposed, in a proof of the 
H-theorem using nonlinear FPEs 1^; as a consequence of such a relation, one obtains 
that the BG entropy is directly connected to the linear FPE, whereas generalizations of the 
BG entropy are associated with nonlinear FPEs [l4|. This reinforces the be 



nonlinear FPEs are intimately related to nonextensive statistical mechanics 
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ief that 



If one considers the nonlinear FPE associated with the nonadditive entropy Sq [iTj in the 



absence of an external force field, one gets 



the porous-medium equation (also known 



2 



as nonlinear heat equation) 



i8i 



dt 



(1) 



which governs the time evolution of the probability distribution P(x, t) for finding a diffusing 
particle in the position x at time t, in a medium characterized by a diffusion constant D; 
the linear diffusion equation comes out as a particular case, by considering q = 1. It should 
be noticed that, within a numerical analysis, as will be one of the main purposes of the 
present work, the parameter q above corresponds to the index associated with the asymptotic 
behavior of the solution of Eq. ([1]), known as a g-Gaussian, related to the nonadditive entropy 
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271 ] ■ These systems consist in assemblies of 



classical rotators, which evolve in tiine, e.g, ., in a plane (XY rotators 
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or in a sphere (Heisenberg rotators 23|, l2J]). For the case of infinite-range ferromagnetic 



interactions, i.e., in the mean-field limit, one may calculate thermodynamic properties an- 
alytically, within the BG canonical ensemble, and in particular, verify the existence of a 
continuous phase transition. The interesting aspect about these systems is that one may 
compute their time evolution through a direct integration of their equations of motion, with- 
out introducing a priori any phenomenological dynamic rules, but just using Newton's law. 
A quite curious behavior has shown up by starting the numerical-integration procedure with 
initial conditions very different from those required in the standard BG equilibrium: the 
system gets trapped in metastable states, before approaching their corresponding terminal 
thermal equilibria. These metastable states are characterized by "kinetic temperatures" 
that are different from the equilibrium ones; besides that, the duration of such states in- 
creases with the number of rotators, A^. Hence, if one considers the thermodynamic limit 
(A^ oo) before the long-time limit, these systems will remain in these metastable states 
and will never reach their terminal equilibrium state, in such a way that the phase space 
will not be equally and completely covered, i.e., these systems are nonergodic. Moreover, 
in such metastable states, the maximum Lyapunov exponent approaches zero, as A^ ^ oo 



20 



23|, contrary to what is expected in a standard BG equilibrium state. For finite values 



of A^, at large enough - but realizable computational times - one approaches a state that 



is presumably the terminal thermal equilibrium state, in the sense that its kinetic tempera- 
ture is in agreement with the one obtained through the BG canonical-ensemble calculations. 
However, although the kinetic temperature of the long-time limit state agrees with the one 
obtained from the canonical-ensemble calculations, other properties may still not coincide 
with those expected in a true equilibrium BG state; as an example, one has a recent analysis 
of the angles described by the infinite-range-interaction XY rotator model, for which their 



distribution in the long-time-limit is well- fitted by a g-Gaussian, with g ^ 1.5 [2l|. This 
suggests that the BG equilibrium state is approached through different steps, one of them 
being the attainment of the equilibrium temperature; a relevant question concerns how long 
will it take for the system to reach completely the final BG equilibrium state. 

In the present work we search for clues on how the approach to equilibrium occurs in 
the above-mentioned Hamiltonian models, by investigating the time evolution of probability 
distributions in a much simpler system, i.e., the porous-medium equation. For that, we in- 
tegrate Eq. ([T]), starting the integration procedure with an initial distribution different from 
its asymptotic solution. In particular, we will consider as initial distribution a g-Gaussian 
characterized by an entropic index qi [qi ^ q), and will follow the time evolution of such a 
distribution towards the final (i.e., asymptotic) g-Gaussian, specified by the index gj, that 
hopefully, g/ = g. Exploring the stability of g-Gaussians in an environment given by Eq. ([1]) 
may help to understand why metastable states that appear in the infinite-range-interaction 
models of rotators remain stable over such long periods. In the following section we discuss 
the exact solutions of the porous-medium equation and their connection to anomalous diffu- 
sion. In Sec. mil we introduce generalized moments, as well as a generalized kurtosis, which 
are more appropriate for dealing with fat-tailed distributions. In Sec. IIVI we analyze the 
time evolution of a probability distribution by following the linear diffusion equation [g = 1 
in Eq. ([1])], provided that the initial distribution is given by a g-Gaussian with gj 7^ 1. In 
Sec. |V]we carry a similar analysis for the porous-medium equation [g 7^ 1 in Eq. ([T])], having 
as initial state a g-Gaussian with g^ ^ q. In Sees. IIVI and IVl we show, by monitoring the time 
evolution of the kurtosis, that the approach to the final g-Gaussian obeys a g-exponential 
function, characterized by a relaxation index q^-d. At this point, it is important to stress 
that in the present work we deal, in principle, with four indexes: (i) the index g defined by 
Eq- dl]); (ii) Qi, associated with the initial g-Gaussian distribution; (iii) qj, associated with 
the final g-Gaussian distribution (for which one expects, g/ = g); (iv) grd, related to the 
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g-exponential of the relaxation towards the asymptotic distribution. However, in the present 
work we have found no evidence, either in our analytical or numerical approaches, of qj ^ q; 
hence, we shall assume from now on that gj = q. Therefore, we will restrict our analysis to 
three indexes, namely, q, qi, and q^^i, as defined above. Finally, in the last section we present 
our main conclusions. 

II. EXACT SOLUTIONS OF THE POROUS-MEDIUM EQUATION 

In this section we discuss briefly the well-known exact solutions for a diffusing particle 
following Eq. ([1]). In order to guarantee the preservation of the normalization for all times t, 
one should impose the probability distribution, together with its first derivative to be zero 
at infinity. 



P{x,t)\ 



x— >itoo 



dP{x,t) 



dx 



, (Vt) . 



(2) 



x— ►±oo 



If one chooses a perfectly localized particle as the initial state, P{x, 0) = S{xo) [S{x) de- 
notes the delta- function] , then, following Refs. j?, Q], one can write the solution of Eq. ([T]), 
satisfying the conditions of Eq. ([2]), in terms of a g-Gaussian, 



P{x,t) = Zgbq{t) eq 



b^Jt)(x-xor 



for g < 3 



(3) 



where e^ = [1 -|- (1 — q)x]^^ (herein, the bracket [C]+ = C, for C > 0, and is zero 
otherwise) represents the g-generalization of the standard exponential function that is re- 
covered in the limit g — ^ 1; its inverse, known as the g-logarithmic function, is given by 



lUq X 



\x 



l-q 



— 1)/(1 — g). This solution presents a compact support for g < 1, and ex- 



hibits power-law tails for g > 1. The time-dependent part of the solution, and the 

normalization constant, Zg, are given, respectively, by 



hq{t)=[2D{2-q){?,-q)Zl~^ t] 



1/(9-3) 



(1^(2 - g) >0) 



(4) 
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where T{x) represents the Gamma-function. 

The diffusion is usually characterized by the time behavior of the second moment of 
the distribution, which is given by < >= scaling as t^l^'^~i\ Hence, q = 1 

yields a linear increase in time, i.e., normal diffusion, whereas Eq. ([T]) leads to anomalous 
diffusion for q ^ 1. Within anomalous diffusion, one may distinguish super-diffusion (g > 1), 
characterized by long-tailed distributions, from sub-diffusion (g < 1), related to compact- 
support distributions. 

If one starts the numerical integration of Eq. ([1]) with P(x, 0) = 5(xo), one follows the 
corresponding g-Gaussian, associated with < >~ t^/^^~i\ However, if one uses as an 
initial distribution a g-Gaussian specified by an index g^ ^ g, the numerical procedure will 
take some time to gradually change from such an initial, to the final g-Gaussian distribution. 
Herein, we will be particularly interested in measuring the time that the system takes to 
approach its final distribution asymptotically. Therefore, different initial states of the system 
are expected to yield different relaxation behavior. In the analysis that follows, the solution 
presented in Eqs. (I3])-(I5]) will be our reference and its kurtosis, to be introduced in the next 
section, will be an important quantity to characterize the relaxation behavior. 

III. GENERALIZED MOMENTS AND KURTOSIS 

The kurtosis is usually defined in terms of the ratio between the fourth, and the square 
of the second moments of a given distribution. However, for the g-Gaussian distribution 
defined in the previous section, one gets divergences in its even moments, in such a way 
that the second moment diverges for g > 5/3, whereas the fourth moment diverges for 
g > 7/5. Therefore, the standard definitions of moments and kurtosis become useless for 
certain ranges of g values. In the present section, we introduce generalized moments, and 
apply them in a definition of a generalized kurtosis. 
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Let us define the generalized n-th moment of a given distribution as, 



< X >, 



J dx x"[P(x,t)]' 

-oo 
oo 

/ dx[P{x,t)Y 



(6) 



wliere n is a positive integer and r > 0. For the g-Gaussian distribution, defined in the 
previous section, one has 



< X >r 
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for < (g - 1) < |r , 

for g = 1 , 
for g < 1 , 

for < (g - 1) < |r 
for g = 1 , 
for g < 1 . 
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(7) 
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Considering the above moments, one may define the following generalized kurtosis. 



< X^ >r 



ro n^-f)(^(^-i)) 

r(i^+f)(r(Tf,+|)) 
rr(^+|)(r(^+|)) 



2 2 \ 



for < (g — 1) < min (g, , ^ 

for g = 1 , 
for g < 1 , 



(9) 



If one uses the property r(a; + 1) = xr{x), it is possible to write this kurtosis in a form that 
covers all three possibilities above. 



«r-,s(g) = 3 



(25-3(g-l))2 



(2r-3(g- l))(2r-5(g- 1)) 



, 2 2 

(g — 1) < min -r, -s 
3 5 



(10) 



In the numerical integration of Eq. ([T]) the initial g-Gaussian, characterized by an entropic 
index gj (gj 7^ g), will evolve in time towards the final g-Gaussian. Obviously, it is desirable 
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FIG. 1: The kurtosis at the beginning of the numerical procedure [Eq. (jlip ]. characterized by 
hnear diffusion in the asymptotic regime {q = 1), Ki.i(gj), is exhibited versus qi (dashed hue). 
Also shown is the kurtosis at the final regime, Kf^q{q) versus q (full line). Notice that the standard 
Gaussian value, ki^i(1) = 3, is recovered from both cases. 

to have Hr,s{(l) always finite when the parameter q varies in the interval — g; unfortunately, 
this is not possible in some cases, and for this purpose, one has to choose the exponents 
r and s conveniently. Herein, we choose these exponents by giving preference to a finite 
kurtosis in the asymptotic limit (t ^ 1). 

In a recent proof of a generalized central-limit theorem, g-Gaussian distributions appear 



as a result of given composition rules 28|, |29|. Inspired by some results of this theorem. 



we will consider the choices s = q and r = f = (g + l)/(3 — g); the choice r = f ensures 
a finite second moment in the asymptotic limit (t ^ 1). Substituting these quantities in 
Eq. ( fTOl) . one gets the following expression for the kurtosis at the beginning of the numerical 
procedure, 

111 -5- 34.(3 -9)1117-35 -54.(3 -9)] • ^ ' 



which is finite, provided that (g^ — 1) < min[(2/5)(g+ l)/(3 — g), (2/3)g]. An interesting 
particular case of Eq. ( fTTl) is the one characterized by linear diffusion in the asymptotic 
regime, i.e., g = 1, yielding /ti,i(gj) = 3 ■ (5 — 3gi)/(7 — 5gj), which leads to a divergence at 
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qi = 7/5, as shown in Fig. [TJ 

The same choices for the exponents r and s yield the kurtosis of Eq. ffTOj) in the asymptotic 
regime, 



3 ^^^2 _ + ll)(5g2 _ 18g + 17) • ^^^^ 

In Fig. [1] we exhibit the kurtosis above versus q, showing that, as expected, it does not 
diverge. 

Therefore, when the g-Gaussian changes between the two entropic indices q^ q, the 

kurtosis evolves in time, changing its behavior between those described in Eqs. flTTl) and 

(fT2l) . respectively; although the kurtosis may be infinite during its time evolution, it will be 
finite in the asymptotic limit. 



IV. q = l AND ARBITRARY INITIAL DISTRIBUTIONS 

In this section we restrict our study to the linear diffusion equation [Eq. ([1]) with q = 1], 
analyzing the time evolution of different initial states given by g-Gaussians [cf. Eq. ([3])], 
characterized by distinct entropic indexes qi. 

As an illustration, we exhibit in Fig [2] the time evolution of P{x, t), starting the numerical 
integration with two typical g-Gaussians, namely, g^ = 3/4 [Fig[2](a)] and g^ = 5/4 [Fig[2]^b)]. 
On the linear scale, one observes no notable difference between the time evolution of these 
probability distributions. Nevertheless, on the log-linear scale (cf. insets of Fig. [2]), one sees 
clearly that the initial distribution relaxes faster to the Gaussian limit in the case g^ = 3/4, 
whereas the fat tails remain stable over a longer period for g^ = 5/4. In particular, in 
this later case, one notices the presence of an inflection point, characteristic of g-Gaussians 
with g > 1, at intermediate times (i.e., in the transient regime), which disappears when the 
distribution approaches the Gaussian limit. 

In what concerns the kurtosis, Ki i(gj), introduced in the previous section, it may be 
calculated analytically, for gj < 7/5. This quantity enables a measure of the time the system 
needs to reach asymptotically the Gaussian distribution; from now on, we will indicate its 
time dependence explicitly, by referring to it as ni^i{qi,t). As we deal here with a linear 
partial differential equation, the calculation of the moments may be carried out exactly by 
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FIG. 2: Time evolution of the probability distribution associated with the linear diffusion equation 
at typical times, t = 0, 0.1, 1, 10 (from top to bottom). The system is initialized with q-Gaussians, 
characterized by initial entropic indexes qi = 0.75 (a) and qi = 1.25 (b). In the insets we exhibit 
the same distributions on a semi-logarithmic scale. 

the use of the Green's function method. The final distribution is a standard Gaussian, 



P{x,t) 



2VM 



exp 



ADt 



(13) 



and the corresponding Green's function is given by 30 1 



10 



Giy,x\t) = ■^^=exp 



ADt 



(14) 



from which one obtains the time- dependent solution of the diffusion equation, 



P{x,t)= / G{y,x\t)P{y,Q)dy 



(15) 



for arbitrary initial distributions P(x, 0). Let us now calculate the n-th moment of the 
time-dependent solution, 



< X >i 



C30 OO OO OO 

j dxx^ j dyG{y,x\t)P{y,Q)= J dyP{y,0) J dx x^G{y,x\t) 

— OO — OO — OO — OO 

OO OO ^ 

= -^1 dyPiy^O) J dxix-yrexp(-^^ . (16) 



Using Eq. (1161) . we obtain the time evolution of the moments of P{x,t) for arbitrary initial 
functions P{x,0), just by calculating standard Gaussian integrals. In order to obtain the 
kurtosis, we calculate the second and fourth moments, respectively. 



< X >i 



OO 

J dyP{y,0) (y2 + 2t) =y2 + 2t 



(17) 



and 



< >i= / dyP{y, 0) {y^ + Uty"^ + Ut^) = y^ + 12yH + 12t^ 



(18) 



where y'^ and y^ denote the standard second and fourth moments of the initial distribution 
P(x, 0). Using these results, the kurtosis becomes. 
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Therefore, the kurtosis' asymptotic value, hmt^oo t) = 3, is approached according to, 

-i,ifet)-3 = ^^^^. (20) 



The equation above yields the kurtosis in terms of the second and fourth moments of initial 
g-Gaussians, provided that qi < 7/5. Notice that i(gj) — 3 is negative, for qi < 1, and 
positive, for 1 < < 7/5. Defining h = bq.{t = 0) = constant, and using Eqs. ([7j) and ([8]), 



fe,<5/3), (21) 



and thus. 



t)-3 = 6 — ^ — = 6 e-ff^'-'''^' , (g, < 7/5). (23) 



The equation above indicates that all initial g-Gaussians, characterized by g^ < 7/5, present 
a kurtosis that will relax to the standard Gaussian (g = 1), following g-exponentials with 
the same relaxation index, grci = 3/2, but different relaxation times, l/[4b'^{5 — 3qi)]. Such a 
dependence of the arguments of these g-exponential functions on the initial entropic index 
qi imply on longer relaxation times for larger values of q^. This is in agreement with the 
results of the numerical calculation exibited in Fig [2], corresponding to the time evolution of 
two distributions initialized as g-Gaussians, with gj = 3/4 and g^ = 5/4, respectively. 

As a test for our numerical algorithm, the result of Eq. fl2^ was reproduced by a numerical 
integration, for typical values of qi, as shown in FigO The integration was carried out using a 



method based on distributed approximating functionals 3l|]- In all cases, we considered the 



distribution of Eq. ([3]) with bg^{t = 0) = 1, as the initial state. The results are represented 
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FIG. 3: Time evolution of the kurtosis obtained by a numerical integration of the linear diffusion 
equation [Eq. ([T]) with g = 1] is compared to the one of the kurtosis Ki,i{qi,t), calculated exactly 
[Eq. ([23j) ]. for typical values of The data is exhibited in (a) linear- linear, (b) log-log and q- 
logarithm-linear [inset of figure (b)] plots. The straight lines in the inset of figure (b) [notice that 
the data for Qi = 0.6 and qi = 0.8 appear essentially superposed] ensure the index gj-ei = 3/2 of the 
relaxation process, and their slopes yield the corresponding relaxation times. 

in different scales, like the linear-linear [Fig. [31(a)] and double logarithm [Fig. El^b)] ones. 
However, an elegant way to show that Ki^i{qi,t) decreases as a g-exponential function, with 
Qvci = 3/2, is by representing the data in terms of the inverse function, i.e., the corresponding 
g-logarithmic function. This is exhibited in the g-logarithm-linear plots in the inset of 
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FIG. 4: Time evolution of the inflection-point position, Xinfi (scaled by Vt), for typical initial q- 
Gaussians, characterized respectively, by qi = 1.2,1.5, and 2.0 (from top to bottom). The small 
fluctuations are just a numerical artifact. 

Fig. El^b), where one observes perfect linear fits, whose slopes give the associated relaxation 
times, l/[4(5 - Sg^)]. 

In a recent analytical work it was shown that the diffusion equation, when initialized with 
a g-Gaussian distribution (1 < < 3), will asymptotically approach its final solution, i.e. 
a Gaussian distribution (g = 



321]. Herein, we present a numerical method to investigate 



how this change occurs for different values of qi within this range. One may see easily that 
the g-Gaussian distributions, defined in Eq. ([3]), present an inflection point for g > 1, when 
represented in a semi- logarithmic plot; the same does not occur for g < 1. Therefore, one 
expects that in the transformation process from an initial g-Gaussian, characterized by an 
entropic index gj > 1, to the asymptotic Gaussian, with g = 1, such an inflection point should 
remain at intermediate times, i.e., in the transient regime, and afterwards, it should approach 
infinity. An example of this effect is shown in the inset of Fig. Et^b), where we exhibit the 
corresponding time evolution of an initial g-Gaussian distribution denoted by g^ = 1.25. 
Hence, the time evolution of this inflection point may provide some additional information 
regarding the process of approach to the Gaussian distribution, and in particular, in the cases 
qi > 7/5, for which the kurtosis of Eq. (l23l) is not defined. Let us herein denote the position of 
the inflection point by Xinfj; we analyze the time evolution of the rescaled quantity, (xinfl/ y/t), 
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in order to measure the time evolution of the inflection point taking off the usual spreading 
effect of the distribution during the diffusion process. We have followed the time evolution 
of (xinfi/v^) for different initial values qi, as shown in Fig. |H Within the time interval 
feasible for computational purposes, we have noticed that {xina/ Vt) always increases in time 
and hopefully diverges, in agreement with the results of Ref. 32]. However, one notices 
that for higher values of qi, such an increase occurs very slowly, and in particular, the case 
qi = 2 suggests that the transformation to the asymptotic Gaussian distribution should take 
place at a very long time. The fact that a single diffusing particle, described in terms of a 
linear equation [Eq. ([T])], may take a very long time to reach its asymptotic-diffusing regime, 
supports the result found in some Hamiltonian systems, described by a set of coupled 



differential equations, for which, given some initial conditions, the final 



never be reached in the thermodynamic limit (A^ — > oo) 
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V. GENERAL CASE: q^l 

In this section we analyze the general case q ^ 1, which corresponds to the nonlinear 
porous-medium equation. The investigation of the solutions of such equation was done 
through a numerical integration of Eq. ([T]), using the same method applied in the previous 
section. Therefore, we followed the time evolution of the kurtosis, which starts from Kf^q{qi, 0) 
[Eq. ( fTTl) ] and will evolve towards its asymptotic limit, Hf,q{q) [Eq- ( p^ ]. In particular, we 
will search for the relaxation law associated with |Kf^q(g) — K,f,gi(li,t)\. 

In Fig. [5] we exhibit the quantity \Kf^g{q) — Hf,q{(li,t)\ for two typical final g-Gaussians, 
namely g = 3/4 and q = 5/4, starting the numerical procedure with different initial q- 
Gaussians. Similarly to what happened in the previous section, our numerical investi- 
gation yields, in both cases, that the kurtosis relaxes to the corresponding final values, 
i^f,q{q) = 2.34797... (g = 3/4) and i^r,q{q) = 3.81717... (g = 5/4), according to g-exponentials, 
whose relaxation index grei depends only on g. The full straight lines in Figs. [5t^a) and ^h) 
correspond to the power-law decays, and which are associated with g-exponentials 
characterized by the relaxation indexes, grci = 7/5 (for g = 3/4) and grci = 5/3 (for g = 5/4), 
respectively. Taking into account the results obtained in the previous section as well, i.e., a 
relaxation following a g-exponential with grci = 3/2 (for g = 1), we propose a general form 
for the relaxation of the kurtosis. 
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FIG. 5: Time evolution of the absolute values of the differences between the kurtosis at their initial 
and final values, i^f,qiQi,t) [Eq. pip ] and Kf^q{q) [Eq. (jl2p ]. respectively, for typical initial and final 
g-Gaussians. (a) q = 3/4 and several initial values of qf, the full straight line corresponds to the 
power-law decay (b) q = 5/A and several initial values of qf, the full straight line corresponds 

to the power-law decay (c) The data of (a) is represented in a log^/s scale, (d) The data of 

(b) is represented in a log5/3 scale. 



l^fM - l^f,,{q^,t)\ = A{q„q) [l + (1 - qrci)b' f{q., q)tY^^' '^"^ , (24) 

where the index gj-ei associated with the g-exponential relaxation process depends only on 
the index q, characteristic of the asymptotic g-Gaussian distribution, and that it appears to 
follow the heuristic relation, 

^-i(^) = ^ ■ (25) 
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FIG. 6: The argument f{qi,q) of the g-exponential characterizing the relaxation of the kurtosis, 
from an initial g-Gaussian (index qi) to an asymptotic g-Gaussian (index q), as defined in Eq. (j24p . 
is exhibited as a function of qi, for different values of g. In the case q = 1, the points (represented 
by squares) were computed numerically, whereas the full line corresponds to the analytical result. 
In the case q = 0.75 {q = 1.25) the points were computed numerically, and the straight dashed 
(dotted) line corresponds to a linear fit. 

The coefficient A{qi, q) that appears in Eq. (1241) should satisfy A{q, q) = 0, in such a way that 
starting the numerical integration procedure with the exact solution of Eq. ([T]), the kurtosis 
should not change in time, i.e., the initial distribution remains stable for all times. The 
argument of the g-exponential, f{qi,q), should recover the particular case q = I, calculated 
analytically, /(g^, 1) = —4(5 — 3qi) [cf. Eq. ( !23l) ]. This argument was estimated numerically 
for typical values of g 7^ 1, as shown in Fig. [6], and our results suggest a simple (essentially 
linear) general form, /(gj, g) = a(g) + b{q)qi. 

VI. CONCLUSIONS 

We have analized, using both analytical and numerical approaches, the stability of g- 
Gaussian distributions as particular solutions of the porous-medium equation. This was 
done by investigating the relaxation towards the final, asymptotic g-Gaussian solution, char- 
acterized by an index g/, when considering as an initial distribution a g-Gaussian, specified 
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by an index g^. By following the time evolution of the kurtosis (defined for qi < 7/5, but 
always finite in the asymptotic limit, t ^ 1), we have found evidence that such a relax- 
ation process follows a g-exponential function, characterized by a relaxation index qreiil)- 
Therefore, in principle, the problem considered may be formulated in terms of four indexes: 
(i) the index q defined by the porous- medium equation; (ii) gj, associated with the initial 
g-Gaussian distribution; (iii) gj, associated with the final g-Gaussian distribution; (iv) q^-d, 
related to the g-exponential of the relaxation towards the asymptotic distribution. Since we 
have found no evidence of qj 7^ g, we assumed that qf = g; this supposition is also supported 
jy a recent analytical approach of the asymptotic behavior of the linear diffusion equation 



32( 1 ■ Accordingly, our study was restricted to three indexes, namely, g, g^, and grei- 

In the definition of this kurtosis, the powers {s and r) that appear in the probability 
distributions of the second and fourth generalized moments were chosen conveniently in order 
to yield a finite kurtosis in the limit t ^ 1; although these choices are arbitrary, we expect 
that other alternatives (e.g., those used in Ref. [33|) should not change the present results 
qualitatively. By using a numerical approach based on the evolution of the inflection point 
that appears in a semi-logarithmic plot of a g-Gaussian with g > 1, we have observed that 
in some cases, an initial infinite- variance distribution (g^ > 5/3) may take a very long time 
to be transformed into a finite- variance one (g < 5/3). In particular, considering the linear 
diffusion equation, we have shown through this method that an infinite-variance distribution 
{Qi > 5/3) evolves very slowly in time towards the asymptotic Gaussian distribution. The 
fact that a single diffusing particle, described in terms of a linear equation, may take a 
very long time to reach its asymptotic-diffusing regime, supports the existence a metastable 
state found in some highly-interacting Hamiltonian systems, described by a set of coupled 
linear differential equations, whose duration diverges in the thermodynamic limit (A^ 00). 
Moreover, for a finite (but sufficiently large) A^, it has been found recently that the angles 
described by the infinite-range-interaction XY rotator model follow a distribution in the 
long-time-limit (i.e., in the limit for which its kinetic temperature coincides with the one 
of the BG canonical ensemble) that is well- fitted by a g-Gaussian, with g ^ 1.5 2l|. This 
suggests that in such Hamiltonian models the BG equilibrium state is approached through 
different steps, one of them being the attainment of the equilibrium temperature; in the 
simpler system considered herein, the approach to the final, asymptotic solution, follows a 
relaxation behavior that may be also very slow in some cases. 
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